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^ ■ The neutralino, the lightest stable supersymmetric particle, is a strong theoretical 

candidate for the missing astronomical "dark matter" . A profusion of such neutrali- 
nos can accumulate near the Sun when they lose energy upon scattering and are 
gravitationally captured. Pair-annihilations of those neutralinos may produce very 
high energy (VHE, above 100 GeV) gamma-rays. 

Milagro is an air shower array which uses the water Cherenkov technique to detect 
extensive air showers and is capable of observing VHE gamma-rays from the direction 
of the Sun with an angular resolution of 0.75°. Analysis of Milagro data with an 
exposure to the Sun of 1165 hours presents the first attempt to detect TeV gamma- 
rays produced by annihilating neutralinos captured by the Solar system and shows 
no statistically significant signal. Resulting limits that can be set on gamma-ray 
flux due to near-Solar neutralino annihilations and on neutralino cross-section are 
presented. 



PACS numbers: 95.35,+d, 14.80.-j 
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INTRODUCTION 



There is very strong evidence that the Universe, and the galaxies in particular, are full of 
non-baryonic "dark matter" (see, for example, [3, 21 , |22[ 0, HH HE 13]). One candidate 
for this dark matter is the neutralino (V) — a weakly interacting massive particle (WIMP) 
predicted by super-symmetric theories 3 3 ■ Experimental tests for this possibility include 



direct searches with extremely sensitive devices which can detect energy deposited by a 
neutralino when it elastically scatters off a nucleus and indirect searches which look for 
products of neutralino-neutralino annihilations. 

The Italian/Chinese collaboration (DAMA) employs an ionization bolometer and has 
reported an observation which they interpret as consistent with the annual modulation 
predicted if WIMPs exist • However, there are possible modulating systematic errors. The 
CDMS experiment detects phonon vibrations of the crystal lattice caused by the neutralino- 
nucleon scattering in the detector volume and has obtained data that appear to exclude 
most of the DAMA-allowed region p| . They reach a spin- independent WIMP-nucleon cross- 
section limit around 2 • 10~ 42 cm 2 in the mass range 20 — 100 GeV. Edelweiss, which utilizes 
both phonon and ionization bolometers, has also released results that significantly cut into 
the DAMA allowed region p. 

Indirect searches look for decay products from neutralino annihilations coming from re- 
gions with enhanced neutralino densities. Searches vary by the regions explored (such as the 
Galactic center or the Sun) and by the decay products being detected. Unlike direct searches, 
interpretation of these experiments requires assumptions about astrophysical parameters 
which give rise to the neutralino annihilation distribution in the region being studied. They 
also depend on cross-sections and branching ratios, both of which are supersymmetry-model 
dependent. 

An increased density of neutralinos may exist in the vicinity of the Galactic center or the 
Sun 23]. This could have arisen in the initial formation of these objects. Also, neutralinos 
entering the Solar system (or the Galaxy) may lose energy via elastic scattering with ordinary 
matter and become gravitationally trapped. Due to the capture and repeated scatterings, 
there would be a near-solar (or Galactic-center) enhancement in the neutralino density. Such 
a local neutralino build-up may provide a detectable flux of annihilation products. Exam- 
ples of indirect searches are those by the Kamiokande and SuperKamiokande underground 
neutrino detectors that have set limits on solar and terrestrial neutralino-induced muon 
fluxes jlE, Such searches are helped by the ability of neutrinos to escape from their 
production region and by the typically large predicted yield, but are hindered by the small 
detection probability. 

Another possible method for detecting dark matter particles is from their annihilation into 
7-rays. The Minimal Supersymmetric extension of the Standard Model (MSSM) predicts 
that the gamma rays emerging from \X ~> 77 an d XX %1 neutralino annihilation modes 
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will give distinct monochromatic signals in the energy range between 100 GeV and 10 TeV, 
depending on the neutralino mass. An additional "continuum" spectrum signal of photons 
will be produced by the decay of secondaries produced in the non-photonic annihilation 
modes. However, past and present high energy gamma-ray experiments, such as EGRET 
and the Whipple atmospheric Cherenkov Telescope, lack the sensitivity to detect annihilation 
7-line fluxes predicted for many allowed supersymmetric models and Milky Way halo profiles. 
The next generation ground-based and satellite gamma-ray experiments, such as VERITAS 
and GLAST, will allow exploration of portions of the MSSM parameter space, assuming 
that the dark matter density is peaked at the galactic center |6[. 

The Milagro 7-ray observatory, which has been taking data since 1999, is sensitive to 
cosmic gamma rays at energies around 1 TeV and is capable of continuously monitoring 
the overhead sky with angular resolution of 0.75°. In this paper, we present the results of 
a search for a TeV gamma-ray signal from the vicinity of the Sun (1-2 solar radii) with 
Milagro. 



II. MILAGRO DETECTOR 

The Milagro Extensive- Air Shower Array is located at 35.88° North latitude and 106.68° 
West longitude in the Jemez Mountains near Los Alamos, New Mexico. At an altitude of 
2630 m above the sea level, its atmospheric overburden is about 750 g/cm 2 . The Milagro 
detector, commissioned in June of 1999, records about 1700 extensive air shower events per 
second and is sensitive to gamma-showers with energies above ~100 GeV. The detector 
consists of a 21 kiloton water-filled pond instrumented with two layers of photo-multiplier 
tubes (PMTs). These PMTs detect Cherenkov light produced by secondary shower particles 
which enter the pond. The top(bottom) layer has 450 (273) PMTs arranged on a 2.8 x 2.8 m 
grid, 1.4 (6) m below the water surface. 

The direction of each shower is determined as the normal of the plane fitted to the PMTs' 
times using an iterative weighted x 2 -method which rejects outliers. The weights for the x 2 - 
fit are prescribed based on the PMT signal strength. The angular resolution of the array is 
estimated to be 0.75°. 

Extensive air showers produced by cosmic rays are the primary source of triggers in the 
experiment. These showers are likely to contain hadrons that reach the ground level and 
produce hadronic cascades in the detector, or muons that penetrate to the bottom layer, and 
will thus illuminate a relatively small number of neighboring PMTs in that layer. Photon 
induced showers, on the other hand, generally will produce rather smooth light intensity 
distributions. Based on this observation, a technique for identification of gamma versus 
hadron initiated showers has been formulated |4j and according to computer simulations can 
correctly select about 90% of hadron initiated showers and about 50% of photon induced 
ones. 

The fluctuations in the shower development, small size of the detector and fluctuations in 
its response make energy determination on the event-by-event basis difficult. The absolute 
energy scale can be determined by examining the displacement of the shadow of the Moon 
due to the Earth's magnetic field |3ol ]. For a more detailed description of the detector itself 
and reconstruction techniques used, see references (H, 0,111]. 
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III. DATA ANALYSIS. 

While it is difficult to tell the difference between cosmic-ray and gamma-ray initiated 
showers, a gamma-ray signal can be detected as an excess of events from the direction of 
the Sun above that expected from the cosmic-ray background. The data analysis procedure 
thus entails determining the average expected background signal N% n , counting the number 
of events from the direction of the Sun N on and determining the statistical significance of 
any excess found. 

For each position of the Sun N% n is found using event rates from the same local region of 
the sky at a time when the Sun is not present using the technique described in jTlJ. This 
method is based on isotropy of the cosmic-ray background and the assumption of short time 
scale detector stability. It allows exclusion of known sources from background estimation 
and correct restoration of the number of excess events to be used for flux measurement. 

The statistic U chosen to test for an excess is: 

N - N b 

U= on °" (1) 
laN on + Nl 



where a is the relative exposure ratio of the " on-source" region to the "off-source" one. 
For detailed discussion on how to deal with varying a see [11]. In the absence of a source, U 
approximately obeys a Gaussian distribution with zero mean and unit variance. 1 If a source 
is present, U will still have an approximately Gaussian distribution with unit dispersion but 
with a shifted mean. 

For the current search we define at the outset the critical value of the statistic U as 
u c = 5 (which approximately corresponds to the level of significance of 2.9 • 10~ 7 ). If no 
excess from the vicinity of the Sun is found, we construct a limit on the source strength as a 
strength which, if present, would have given us a detectable signal (u > 5) with 97.7% (2a) 
probability [9]. In addition, following the more standard procedure assuming that the 
source exists, a 90% one-sided confidence interval on its strength is also calculated. 

Given the detector response to particles of different types and assumed source features, 
it is possible to predict the number of events N k to be observed due to the source. 



N k = [ F(k,E)S(k,Q)T(Q)\[ _ A k (E, 0, 0) dO 
J Ueen 



dE dQ (2) 

'sen J 

where F(k, E) is the number of particles of type k with energy E emitted by the source 
per unit area per unit time, S(k,Q) describes the known geometrical shape of the source and 
is assumed to be energy independent, T(Q) is the time during which the source is located in 
local direction and is the direction of the particle arrival as output by the reconstruction 
algorithm. The function A k (E, 0,0) is known as the "effective area" of the detector 2 and 
can be integrated over prior to data processing for selected configuration of Cl. Therefore, 
after integration over 0, the number of events due to particles of type k to be observed from 
a region Cl is: 



1 The conditions of applicability of the Gaussian approximation [ll| are satisfied for the presented analy- 
sis 0. 

2 The Milagro detector response function Jq e q Ak{E, 0, O) d(~)J was obtained by simulating the operation 
of Milagro detector with CORSIKA and GEANT simulation packages (see |4|). 



5 



N k = J F(k,E)A T k (E)dE (3) 

The integrated effective area A\ (E) is obtained during the data processing from Monte 
Carlo generated effective area A k (E, 6, 6). By counting the number of excess events in an 
observation bin iV 7 = N on — N% n , it is thus possible to deduce some properties of the source 
function F(^, E). 

IV. EFFECT OF THE SUN SHADOW 

The gamma ray signal from neutralino annihilations near the Sun should appear as an 
excess number of events from the direction of the Sun over the expected cosmic-ray back- 
ground. The interpretation of any observed signal, however, is not an easy problem. Largely, 
this is due to the fact that the cosmic-ray background is not expected to be uniform; the 
Sun absorbs the cosmic rays impinging on it and forms a cosmic-ray shadow. The situa- 
tion is complicated by the magnetic fields of the Earth and the Sun. Due to bending of 
charged-particle trajectories in magnetic fields, the Sun's shadow in the TeV range of par- 
ticle energies will be smeared and shifted from the geometrical position of the Sun. On the 
other hand, in the presence of strong Solar magnetic fields, lower energy particles cannot 
reach the surface of the Sun and are reflected from it. Such particles are not removed from 
the interplanetary medium and, since the cosmic rays are isotropically incident on the Sun, 
may not even form a cosmic-ray shadow. In addition, the Solar magnetic field varies with 
time, which will smear the shadow in a long-exposure observation. Therefore, it is difficult 
to ascertain the exact shape of the cosmic-ray shadow at the Sun's position and deduce an 
excess above it. 

Nevertheless, the effect of the Earth's magnetic field and the Solar wind can be studied 
by observing the shadow of the Moon during the solar day. The effect of the geomagnetic 
field on the shadows of the Sun and the Moon should be very similar since the Sun and 
the Moon cover similar size regions on the celestial sphere and traverse similar paths on the 
local sky in one year of observation. In addition, the Earth's magnetic field at the Moon 
distance is already so small that any additional deflection beyond the Moon distance by this 
field of particles originating from the Sun can be neglected. 

A deficit of events due to such a shadow could be filled by a signal from neutralinos 
which would remain undetected when the excess is searched for. Since the shadow depth 
is unknown, when performing a search for evidence of a positive excess we are forced to 
assume the absence of the shadow. On the other hand, for setting a conservative limit on 
the photon flux we use the maximum possible depth of the shadow. 

The deficit of events from the direction of the Sun can not be greater than that produced 
by the Moon because as mentioned above the additional Solar magnetic field can only 
decrease its shadow relative to the Moon. To be conservative in setting the gamma-ray 
flux limit, the strongest event deficit produced by the Moon in 5° radius from its position, 
corrected for relative Sun/Moon exposure, is used as a correction for the possible presence 
of a shadow of the Sun. 
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V. RESULTS 

A. Outcome of the observation 

The data used in this work was chosen to satisfy the following conditions: online recon- 
struction between the 19th of July 2000 and the 10th of September 2001, the number of 
PMTs required for a shower to trigger the detector greater than 60, the number of PMTs 
used in the angular reconstruction greater than 20, zenith angles smaller than 45 degrees, 
and passing the gamma/hadron separation cut 0]. The start and end dates correspond re- 
spectively to introduction of the hadron separation parameter into the online reconstruction 
code and detector turn-off for scheduled repairs. Several data runs were removed from the 
dataset which included calibration runs and the data recorded when there were online DAQ 
problems. 

For the solar region analysis, ±5° regions around the Moon and the Crab nebula were 
vetoed from the data set as they present known sources of anisotropy to the cosmic-ray 
background. 

Photons produced in the Sun will be absorbed, whereas the distribution of neutralino 
annihilations outside is a rapidly falling function of distance from the Sun. Therefore, we 
believe that the gamma-ray signal is produced mostly between 1 and 2 solar radii and treat 
the gamma-ray source as circle of 0.5° radius. It has been shown that the optimal bin 
size is a slow function of the source size and for estimated 0.75° angular resolution of the 
detector, the optimal "on-source" bin is a circular one with the radius of 1.26° centered on 
the Sun. 

Overall, 1164.7 hours of exposure to the Sun is obtained in the data set. The total number 
of events observed in the "on-source" bin is N on = 137211 while N% n = 137728 events is 
expected based on the "off-source" exposure, leading to the value of the test statistic U of 
—1.35 (see figure P). Therefore, the null hypothesis of the absence of gamma-ray emission 
from the Sun can not be rejected and a limit on the possible 7-ray flux from the solar region 
is obtained. 

Overall, 423.5 hours of exposure to the Moon during the day time is obtained in this 
data set after ±5° regions around the Sun and the Crab nebula were vetoed. The largest 
deficit observed in the sky map centered on the geometrical position of the Moon is —3.3 
corresponding to —766 events (see figure Q). The exposure on the Sun is 2.75 times greater 
than that on the Moon during solar day, leading to estimated maximal deficit in the Sun's 
direction as: -766 ■ 2.75 = -2107. 

Hence, using the criteria outlined in section 11111 we conclude that the number of excess 
events due to a possible Solar source of VHE photons does not exceed N exdude = 4791. It 
is this number which is used to construct the exclusion region on gamma-ray flux strength 
from the Solar region (equation [SJ) . 

The 90% one-sided confidence interval on the average number of excess events iV con /jrf ej ™ 
assuming that gamma-ray emission from the Sun exists constructed in a standard way [8J] 
is A? confidence — 2081. It is this number which would allow estimation of the gamma-ray 
flux strength from the Solar region if it were known that such gamma-ray emission exists 
(equation EJ). 



(a) 



(b) 
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FIG. 1: The maps of statistical significance of the number of excess events from the regions of the 
sky around the daytime Moon (1(a)) and the S un (|l(b) ) and the corresponding source exposure 
as function of zenith angle in hours per degree ( l(c)|l(d) ). The color code is the value of U (see 
equation lfl|)). The Moon shadow is much less significant than that reported in because of 
7-selection and other cuts used in the present analysis. The maps are made using azimuthal equal- 
area projection in polar case centered on the corresponding Celestial object. Contours represent 
loci of distance 1.26°, 2.0° and 5.0° from the Celestial object. 
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FIG. 2: Integrated effective area A\ as a function of photon (solid) and proton (dashed) energy 
(see equation (J3J) for the circular 1.26° bin centered on the Sun for the trigger condition and cuts 
used in this analysis. 



B. Limit on the photon flux due to neutralino annihilations 

Because the two close direct-production spectral lines can not be resolved by the Milagro 
detector, the differential photon flux due to neutralino annihilations is assumed to have the 
form (see 

dF(E) F f^)~ 3/ V 7 - 8 ^ 

= F 5 5(E y - mJ) + — • ^ ; (4) 

where F$ is the integral flux due to a 5-function-like photon annihilation channel and F c 
is the integral flux of photons with energies greater than 0.01 • m x due to continuum photon 
spectrum annihilation channel of neutralinos with mass m x . Here, the normalization of the 
continuum photon spectrum has been written out explicitly. 

Using this expression for the flux with the integrated effective area of the detector, we find 
a relationship between the number of observed events iV 7 and the integral flux parameters 
Fs and F c in the form: 

FgA + F C S = iV 7 (5) 
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m =0.2 TeV 



m =0.5 TeV 




m v =2.0 TeV 



m v =5.0 TeV 




F a , cm s 



F t , cm s 



FIG. 3: The values of (F$,F C ) below the lines are allowed for corresponding neutralino masses. 



m x (TeV) 


A (cm?s) 


S (cm?s) 


0.1 


1.055 x 10 11 


0.000 


0.2 


8.772 x 10 11 


4.969 x 10 7 


0.5 


6.070 x 10 12 


2.634 x 10 9 


1.0 


3.389 x 10 13 


2.127 x 10 10 


2.0 


1.600 x 10 14 


1.280 x 10 11 


5.0 


5.942 x 10 14 


1.208 x 10 12 


10.0 


1.608 x 10 15 


5.575 x 10 12 


20.0 


3.684 x 10 15 


2.136 x 10 13 


50.0 


8.030 x 10 15 


1.035 x 10 14 



TABLE I: Coefficients of the flux limit calculation (see equation Q) 



The values of coefficients A and £ (see table |TJ) are obtained by substituting the expression 
for flux (j3J) and the integrated effective area (figure EJ) into the formula for the number of 
expected events (J3J) - 3 Note, that depending on m x the integral in (j3J) may extend below the 



0.01m x normalization. 



The figure El presents the curves demarcating the allowed and excluded regions in the 
photon flux parameter space (Fg, F c ) corresponding to the significance 2.9 • 10~ 7 and the 
power of the test of 97.7% obtained by setting iV 7 = N exdude which have a form of straight 
lines in the (Fg,F c ) plane. Because there is only one equation (eq. (JHJ)) constraining two 
parameters, table |H] provides the most conservative limits on each of F$ and F c when the 



Should an emission model other than the one presented by equation Q become a theoretical preference, 
figure |21 can be used to recompute the model's parameters. 
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m x (TeV) 


Fs < (cm~ 2 s _1 ) 


F c < (ernes' 1 ) 


0.1 


4.54 x 10~ 8 


— 


0.2 


5.46 x 10~ 9 


9.64 x 10~ 5 


0.5 


7.89 x 10~ 10 


1.82 x 10~ 6 


1.0 


1.41 x 10~ 10 


2.25 x 10~ 7 


2.0 


2.99 x 10~ n 


3.74 x 10" 8 


5.0 


8.06 x 10~ 12 


3.97 x 10~ 9 


10.0 


2.98 x 10" 12 


8.59 x 10" 10 


20.0 


1.30 x 10~ 12 


2.24 x 10" 10 


50.0 


5.97 x 10~ 13 


4.63 x 10~ n 



TABLE II: The limit on the flux parameters (Fs,F c ) corresponding to the significance 2.9 • 10 7 
and the power of the test of 97.7%. For the 90% one-sided confidence interval multiply these limits 
by 0.4344. 



m x , {TeV) 


I(m x )x 10-«cm*0.3GeV/cm> ( X) 


0.1 


1.65 x 10 18 


0.2 


4.17 x 10 17 


0.5 


6.72 x 10 16 


1.0 


1.68 x 10 16 


2.0 


4.22 x 10 15 


5.0 


6.72 x 10 14 


10.0 


1.69 x 10 14 


20.0 


4.22 x 10 13 


50.0 


6.75 x 10 12 



TABLE III: Capture rate I(m x ) as a function of neutralino mass m x obtained in a three-dimensional 
calculation The normalizations of a vx and po are typically used numbers. 

contribution from the other is set to zero. 

For the depiction of the one-sided 90% confidence interval in the (Fs, F c ) plane, both axes 
in the figure H should be rescaled by N exdude /N confidence = 0.4344. 

Figure El is the derived limit on the values of the parameters of the gamma-ray emission 
model (equation (@J)) due to near solar WIMP annihilation and is independent of the models 
of their distribution in the Milky Way galaxy and the Solar system. 

C. Neutralino limits 

The interpretation of the constructed limit on the gamma-ray flux is highly model depen- 
dent. It is based, for instance, on assumptions regarding the shape of the velocity distribution 
of the dark matter in the galactic halo and its density profile in the Solar System. Therefore, 
several assumptions are made to construct limits on physically interesting quantities. 

We assume Maxwellian velocity distribution of the neutralinos in the Solar system with 
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m x (TeV) 






(J/ 1 ) i -ii cm -2 o.ZGeV/cm :ib jf° ut < 


(J /I) 10 -4:i cm 2 0.3GeV/cm^lf° ut < 


0.1 


— 


155 


0.2 


1.30 x 10 6 


73.7 


0.5 


1.52 x 10 5 


66.1 


1.0 


7.54 x 10 4 


47.3 


2.0 


4.99 x 10 4 


39.9 


5.0 


3.32 x 10 4 


67.5 


10.0 


2.86 x 10 4 


99.2 


20.0 


2.99 x 10 4 


173 


50.0 


3.86 x 10 4 


497 



TABLE IV: The limits on J(m x ) / T(m x )a PX p b^f out and J(m x )/ T(m x )a px pob^f out . 
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FIG. 4: The values of yPocr px b^f ou tJ(m x )/I(m x ), Po^ PX b^f ou tJ(m x )/I(m x )J below the lines are 
allowed based the constructed limit for corresponding neutralino masses. As one progresses from 
low to high m x , the detector effective area goes up at the same time the flux of incoming neutralinos 
goes down as l/m x for fixed dark matter density pq. In addition, the capture probability is also 
decreasing as l/m x in an elastic scattering on a fixed mass target. This explains qualitative 
behavior of these lines. 
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FIG. 5: Upper boundary of a 90% one-sided confidence interval on neutralino-nucleon cross-sections 
obtained from different direct-search experiments. The closed contour is the allowed region at 3a 
confidence level from the DAMA experiment. The plot is adopted from ^jj. The two Milagro curves 
plotted assume an equilibrium situation where the annihilation rate J(m x ) equals the capture rate 
I(m x ), and consider the direct and nondirect annihilation modes for preset values of 6 7 

EH. 

The Milagro limit is obtained by dividing the plotted values by the fraction of annihilations that 
occur outside the Sun, f outl for which calculations vary from 10 _1 to 10 



mean velocity Vq = 220 km/s and width 2V 2 . The photon flux of equation (@J) at the Earth 
due to neutralino annihilations can be computed as: 

dF x( E i) = fout(m x ) ■ f e ■ b 1 {E 1 , m x ) ■ ^-yir dE i ( 6 ) 

where b 1 (E 11 m x ) is differential photon yield per neutralino for producing a photon with 
energy E 1 in neutralino-neutralino annihilation, fouti^x) * s ^ ne fraction of neutralinos an- 
nihilating outside the Sun, f e is the fraction of produced photons which escape from the 
Sun and is of the order of 1/2, J{m x ) is the total neutralino annihilation rate in the Solar 
system and = 1.5 • 10 11 m is the mean Sun-Earth distance. 

Given the functional form of the flux from equation (JU), the photon yield b^(E^, m x ) is: 

tft \ ~ 3/2 e -7-8E/m x 

b,(E y ,m x ) = fymJSfa - m x ) + • g^Z^ (7) 

where b^ {b^) is the number of photons produced per annihilation directly (indirectly). 

One can also make an assumption that an equilibrium situation has been reached and 
that the annihilation rate J(m x ) and the capture rate I{m x ) are identical. (We use this in 
figure El below.) 

A 3-D calculation has been performed 0] to determine the rate I(m x ) of WIMP capture 
by the Sun as a function of the neutralino mass. For given local galactic dark matter density 
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p , a structure-less 4 x~P elastic cross-section a px determines how often a neutralino passing 
through the Sun scatters and loses enough energy to get gravitationally captured. The 
results of this calculation are presented in table IIII1 Since I(m x ) is proportional to Po<J px , it 
will be normalized to the capture rate computed at p = 0.3 GeV/cm 3 and a px = 10~ 43 cm 2 . 

A limit on cr px o 7 would provide constraints on parameters of Super Symmetric models. 
Using the formulae ()5I6I7|) . however, one obtains neutralino-mass-dependent limits on the 
product of Po(Tp X bsyf ou t(J/I) as presented in table II V I and figure |U While the value of the 
local dark matter density po is known relatively well, there are substantial disagreements on 
the fraction of neutralino annihilations near the Sun f out . 

The problem of WIMP capture on bound near-solar orbits was considered in [k| 
To our knowledge, the first one-dimensional computer simulation of the distribution of the 
annihilation points near the Sun was treated in |28|. There, it was assumed that the Solar 
system consists of a uniform density Sun only and that the capture of a particle happens in 
its first scattering inside the Sun with an additional assumption that the Solar system has 
reached a dynamic equilibrium and that the capture rate I{m x ) is equal to the annihilation 
one J{m x ). 28] concluded that the fraction of all annihilations happening outside the Sun 
is f out ~ 10~ 5 — 1CT 7 . Under essentially the same assumptions, a simulation done by [I?} 
provided a drastically different prediction of f out ~ 10~ 14 — 10~ 16 . However, a 3-D computer 
simulation of the neutralino annihilation distribution 10] provides an estimate f out ~ 10 _1 . 

Therefore, in figure 03 we illustrate 90% one-sided confidence intervals for two different 
simplified cases: first b s = 0.001 and & 7 = 0, and second 6* = and & 7 = 1. In both cases we 
take po = 0.3 GeV/cm 3 , J(m x ) = I(m x ), and f out = 1. It should be noted that while b\ 



take po = 0.6 uev/cm , J[m v ) = l{m x ), ana j out = i. it snouia De noted tnat wnne ox is 
always less than unity 29] , b c can be as high as 10 for some Super Symmetric models 12 • 



VI. CONCLUSION 

This work presents the first attempt to detect TeV photons produced by annihilating 
neutralinos captured into the Solar system. Analysis of the Milagro data set collected during 
2000-2001 shows no evidence for a gamma-ray signal due to such a process. The limit on 
the possible gamma-ray flux due to such a process with significance 2.9 • 10~ 7 and the power 
97.7% has been set (see table ITD and figure OJ). Even in the absence of a clear signal the 
constructed exclusion limit may constrain the values of free parameters of supersymmetric 
models (see table|TV]and figureEJ). In addition, a standard 90% one-sided confidence interval 
on the magnitude of the photon flux due to near-Solar neutralino annihilations has been 
constructed. 

The interpretation of the constructed limit on the gamma-ray flux is highly model de- 
pendent. Conversion of the flux measurement to a cross section limit requires knowledge of 
the annihilation rate, J(m x ), the fraction of annihilations outside the Sun, f out , and photon 
yield per annihilation, b T These, in turn, depend on parameters of astrophysical and Super 
Symmetric models. As mentioned in the paper, some of these can be estimated using simple 
assumptions. For example, one may assume an equilibrium situation when capture rate 
is equal to the annihilation one. However, absence of a reliable estimate on the fraction 
of annihilations happening outside the Sun makes it hard to interpret the limits in terms 
of the theoretically interesting o 7 o> x . Therefore, the results are presented with all these 



4 Structure-less scattering is the one with no restrictions other than energy and momentum conservation. 
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parameters written out explicitly. 

Continuous improvements in reconstruction algorithms, detector modifications and longer 
observation times will lead to a better upper limit. One of the factors which led to a 
deterioration of the constructed upper limit is the inability to compensate for presence of 
the Solar cosmic-ray shadow due to the intricate structure of the Solar magnetic fields. Once 
the cosmic-ray shadow of the Sun is understood quantitatively, it may be possible to improve 
upon the limit. 
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